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O ; Abstract 

Any exact eigenstate with a definite momentum of a many-body Hamil- 
tonian can be written as an integral over a symmetry-broken function $. For 
Jfl^ ' two particles, we solve the problem exactly for all energy levels and any inter- 

o 

\^ ' particle interaction. Especially for the ground-state, $ is given by the simple 

o' 

■^ ' Hartree-Fock/Hartree ansatz for fermions/bosons. Implications for several 

I' 






'i 



and many particles as well as a numerical example are provided. 
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The homogeneous gas of interacting fermions and bosons is a fundamental concept in the 
physics of many-particle systems, see e.g., Refs. [1,2]. Expect for specific cases, such a many- 
or several-body problem cannot be solved exactly, i.e. the wavefunctions and energies of the 
ground and excited states are not known. In fact, for more than a few particles it becomes 
already impossible to numerically compute the exact ground-state energy and wavefunction. 
Consequently, approximations are a must. 

In this Letter, we introduce an exact continuous configuration-interaction ansatz for the 
many-body wavefunction ^ of interacting particles in a volume with periodic boundary con- 
ditions. Examples for realizations of this case are a ring in one dimension (ID), a torus, a 
long, thin pipe (tube) or a sphere in 2D, and the text-book example of a "big box" in 3D. 
Specifically, we employ a many-body function $ as a basis function for a continuous expan- 
sion of \1'. The shape of $ is to be optimized by employing the variational principle. For two 
particles, we solve the problem exactly for all energy levels and any inter-particle interaction. 
Especially for the exact ground-state, $ is given by the simple Hartree-Fock/Hartree ansatz 
for fermions/bosons. For more particles as will be explained below, our ansatz with any spe- 
cific choice for $ is better than solving the Schrodinger equation by optimizing \1' = <l> itself. 
For instance, taking for $ a Hartree-Fock/Hartree ansatz, the resulting equations would 
lead to lower energies than the corresponding Hartree-Fock/Hartree equations. These prop- 
erties make our ansatz a particularly attractive approximation for the few as well as the 
many-body problem. 

Consider the generic many-body Hamiltonian describing A^ fermions or bosons in a 3D 
box of volume V = L^LyLz'. 

N ff. N 

H{^ur„...) = -Y.Trr2+ E f^(r.-r,), (1) 

Here, Tj = {xi,yi,Zi) (we use Cartesian coordinates) is the coordinate the i-th particle 
and f/(rj — r^) describes a (/enera/ pairwise interaction between the i-th and j-th particles. 
Periodic boundary conditions in 3D are assumed. 

The basic question we would like to address here is how to approach the solutions of 



Eq. (1). Since continuous translational symmetry exists in each dimension, we construct a 
continuous self- consistent configuration-interaction (CCI) ansatz for tlie many-body wave- 
function. Namely, we define tlie CCI wavefunction as 

^(ri, r2, . . .) = / dro C(ro) $(ri - Fq, r2 - Tq, . . .), (2) 

wliere both tlie weight C(ro) and the many-body function $(ri, r2, . . .) are to be determined 
from the requirement that the variation of the energy per particle £{N) = j^ ^^L^^ van- 
ishes, tq is a vector in 3D. Demanding stationarity of £{N) with respect to C(ro), we proved 
that any variational solution ^ must be an eigenfunction of the momentum, say with the 
momentum Pq, where Pq = 27r(^,^,^j and the n^ take some specific integer values. 
This can most generally be fulfilled by fixing C(ro) = e^*"^"'"" and, thereby, allowing $ to 
be a symmetry-broken function (we use this term to describe a function which is not an 
eigenfunction of the momentum). The decomposition of \1' in terms of a symmetry-broken 
$ provides many more degrees of freedom for optimization (this will be evident below) , and 
will lead to lower energies and better approximations for the wavefunctions. In fact, any 
exact solution of Eq. (1) with the momentum Pq can be decomposed into a symmetry-broken 
function $ (such a decomposition is not unique). This finding is central to this work. The 
proof of this property is straightforward to make. The integration over fq eliminates from 
$(ri — To, r2 — To, . . .) any component which is not an eigenfunction of the momentum with 
the eigenvalue Pq. Of course, ^ should possess the other symmetries of the many-body 
problem (1). Specifically, it should be antisymmetric or symmetric under permutations of 
any two particles for fermions or bosons, respectively, and be an eigenfunction of the total 
spin. 

For two particles we can explicitly determine $ which reproduces the exact solution of the 
Schrodinger equation for all energy levels and any interaction U. Specifically, we will prove 
below that taking for $ in Eq. (2) an Hartree-Fock/Hartree ansatz for two fermions/bosons 
is exact for all eigenstates of (1) with Pq = 0, among which is the ground state. The case of 
Pq 7^ is also exact with another "simple" form of $ and will be presented thereafter. We 



refer in the following to the spatial degrees of freedom alone also for two fermions since the 
spin and spatial parts of the wavefunction are separable in this case. A two-particle (spa- 
tial) wavefunction symmetric under permutations of particles' indices, ^+(ri,r2), describes 
bosons and also fermions in the singlet state, whereas an anti-symmetric wavefunction, 
^_(ri,r2), describes fermions in the triplet state. Accordingly, the general solution (with 
Po = 0) takes on a simple form 

^+(ri,r2) = - Y. ^kCos[k-(ri-r2)], 

^_(ri,r2)= E ^kSin[k-(ri-r2)] (3) 

k,A;a;>0 



where k = 27r (^, ^, ^j is the momentum and the rii take integer values. The coefficients 
Ak are real and depend on the nature of the two-body interaction in Eq. (1) and, of course, 
on the energy level and particle statistics. Next, we introduce the following two orthogonal 
orbitals: 

0i(r) = ^= J2 \A4kCos(k-r), 

VV/ k,fe>0 

02(r) = ^= E \/^sin(k-r). (4) 

After some algebra, it is straightforward to show that with these orbitals the CCI ansatz 
(2) is exact, namely that the following decompositions with Hartree and Hartree-Fock wave- 
functions exist: 

^+(ri, r2) = / dro(f)i{ri - ro)0i(r2 - tq), 

^_(ri,r2)= / (irodet |0i(ri -ro)02(r2 -ro)|. (5) 

Let us discuss the results just obtained for Pq = 0. The CCI orbitals (4) are very 
unique. They do not possess the translational symmetry of the problem, i.e., they are 
symmetry-broken functions, even for very weakly-interacting particles, for which the Hartree 
and Hartree-Fock orbitals are symmetry preserving [1,2]. They are, however, eigenfunctions 
of the space inversion operator, r -^ — r. Another interesting property is that 0i(r) and 02(1") 



are, in general, complex quantities. Nonetheless, the resulting exact wavefunction is real and 
possesses the correct translational symmetry. In other words, these additional degrees of 
freedom the CCI ansatz allows for its orbitals make it exact with an Hartree-Fock/Hartree 
ansatz for $(ri,r2), no matter what the inter-particle interaction is. Notwithstanding, 
applying these mean-fields themselves for two particles only is an inaccurate approximation, 
to say least. 

Next, we turn to lay out the case of Pq 7^ 0. Here, the wavefunction may be written as 

k 

(Ak are complex coefficients). Introducing the two general orbitals 
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02(r) = 47E/4;^e^'^^ (7) 

yv k 

the CCI ansatz is again exact with the following decomposition 

^±(ri, rs) = j rfroe+'P""« [0i(ri - ro)02(r2 - tq) ± H^i - ro)0i(r2 - tq)] . (8) 

This concludes our proof that the CCI ansatz with a "simple" choice for $(ri, r2) is exact for 
all two-particle states, for any inter-particle interaction and for fermion and boson statistics. 
Having solved the problem for A^ = 2, we turn to the cases of several and many particles. 
Here we do not expect that an Hartree-Fock/Hartree ansatz for $ will solve the problem 
exactly. That there is such a $ for any solution of (1) was proved above. So, what can be 
said about an approximate $? For any ansatz for $ the variational principle ensures that the 
ground-state energy obtained by optimizing $ within the CCI ansatz is at least that obtained 
by utilizing $ itself as the ansatz for \I'. In practice, the CCI energy will be lower (as found 
above for N = 2), since the CCI ansatz has many more degrees of freedom for optimizing 
$ and minimizing the energy, also see example below. For instance, taking the Hartree- 
Fock/Hartree ansatz for $ the CCI ansatz will provide lower energies than the Hartree- 
Fock/Hartree equations. In conjunction with that, the CCI orbitals will be symmetry-broken, 



even for very weakly-interacting particles, a situation found above for the exact solution of 
the two-particle problem, see Eq. (4). Contrast this with the Hartree-Fock/Hartree orbitals 
which are symmetry preserving for not too strongly interacting particles [1,2]. All the above 
make the CCI ansatz particularly attractive to apply both for many- as well as for several- 
particle systems already with a Hartree-Fock/Hartree ansatz for $. Of course, other ansatzs 
for $ can be plugged into Eq. (2), if needed. 

In the following, we would like to implement the CCI ansatz (2) in its simplest form. 
We concentrate on bosons, set $ to be an Hartree product, <l>(ri,r2, • • •) = HiLi 0(i"i)) P^t 
C(ro) = 1 for the ground state, and derive the CCI equation for the (single-particle) orbital 
0(r). For the inter-particle interaction in (1) we take the contact interaction f/(rj — r^) = 
Uo6{ri — Tj). This interaction is widely used for the dilute Bose gas, see, e.g., Refs. [3,4] and 
references therein, where Uq is proportional to the bosons' s-wave scattering length. 

The final result for e{N) takes on the following appealing form: 



with the energy density (T 



5(iV;ro) 






5(ro) ■ (0(r) T 0(r - ro)) + ^^ '- (0^(r) | 0^(r - ro; 



UoiN-l) ,,^ ,,^2, 



(96) 



The Gross-Pitaevskii (GP) energy [3,4] can be seen as a private case of the CCI energy (9) 
if we were to evaluate the latter for Tq = only. 

By minimizing e{N) in (9) with respect to the orbital 0(r) we obtain an equation for 
this orbital. The final expression reads 

Jdr,S''-\ro) [Siro) ■ f + (N - l)f/o0*(r)0(r - ro)] 0(r-ro) =|rfro5^-^(ro)/i(ro)0(r-ro), 

(10a) 
with the chemical potential density 

/i(ro) = S\ro)Ne{N) - {N - l)e{N - 1; ro). (106) 

Eq. (10), the CCI equation as we shall call it, is an integro-differential equation for the 
orbital 0(r) in the case of A^ interacting bosons, which has to be solved self-consistently. 
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The non-trivial dependence of the CCI equation on the number of particles A^ is via powers 
of the overlap 5'(ro). It comes from the CCI ansatz which allows the function $ and, hence, 
the orbital 0(r) to be translated throughout space. We will see the effect of this dependence 
in a numerical example below. 

Having derived the CCI equation Eq. (10) we would like to solve it for a specific example, 
namely for the attractive Bose gas in ID. This system has attracted a renewed interest 
recently, see, e.g., Refs. [6-8] and references therein. As mentioned in the introduction, in 
ID we may equivalently speak of bosons on a ring. In this case, it is convenient to re-scale 
the coordinate x and define the angle ^p = j^x G [— tt, tt). The effective particle-particle 
interaction becomes then Uo5{(pi — ^2), Uq = ^^^ and T = —-§-^- The properties of the 
system are described by a single dimensionless parameter 7 = °^^~ ' [7] , which is negative 
in the attractive case. 

Let us briefly describe the properties of the GP solution, as far as they are needed here. 
Here, for weakly-interacting bosons, i.e. for I7I < 0.5, the ground-state orbital is a constant 
equal to l/v^27r [6,7], also see Fig. lA. For a stronger interaction, i.e., for I7I > 0.5, the 
GP solution lowest in energy becomes localized and, hence, symmetry-broken [6,7], also see 
Fig. IB. In both regimes, the shape of the GP orbital does not depend on A^ for a flxed value 
of I7I [7]. In Ref. [8] it has been shown that when this fixed, symmetry-broken GP orbital 
is plugged into a CCI ansatz the energy obtained is much lower than the GP one, and that 
"with increasing number of particles and/or strength of inter-particle interaction, is even 
lower than that accessed by tractable diagonalization of the many-body Hamiltonian." In 
this context, Ref. [8] is a special case of this work; There, a symmetry-broken orbital was 
"supplied" by the mean-field solution itself. The breakthrough which we present here is the 
successful optimization of the CCI orbital (j) by solving Eq. (10) for bosons on a ring. To 
this end, we employed the discrete variable representation (DVR) method [5] and solved 
Eq. (10) iteratively (self-consistently) till convergence. The numerical implementation was 
straightforward. Note that, although the wavefunction strongly depends on A^, the numerical 
effort does not. 
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Analyzing the results, we start with the case of A^ = 2 bosons. The Hartree (GP) ansatz 
itself is not a good approximation for such a small number of bosons. Yet, it becomes exact 
inside the CCI ansatz (2), see proof above and Eq. (5). For this, we examine the unique 
shape of the CCI orbital 0. For attractive bosons, this orbital for the ground state is a 
real function (for repulsive bosons this orbital is a complex function). The shape of for 
A^ = 2 is presented in Fig. 1 for the two representative values I7I = 0.2 and I7I = 1.0. The 
corresponding GP orbitals are depicted for comparison. As anticipated on the basis of the 
exact analysis, is symmetry-broken also for values of I7I for which the GP orbital is not. 
For any value of 7 the CCI orbital is totally different from the GP one: it is narrower and 
on top of that exhibits a central spike. Interestingly, the structure of below and above 
I7I = 0.5 is different. This shows that the "phase transition" occurring at I7I = 0.5 manifests 
itself not only within the mean-field (GP) solution [6,7]. Rather, it also appears in the exact 
solution represented by the CCI orbital depicted in Fig. 1. 

Finally, we set on to study the shape of the CCI orbital for N > 2, i.e, the evolution 
of with increasing particle number as obtained by solving Eq. (10) in ID. Here we do 
not expect that employing the Hartree ansatz for $ will provide the exact energy. It will, 
however, provide a very good estimate for it, lower, naturally, from the already low energy 
obtained with the fixed-unoptimized, symmetry-broken GP orbital [8]. In Fig. 2 we present 
for A^ = 5,25,100,1000 and 10000 and I7I = 0.2 the CCI orbital, which is, as expected, 
symmetry-broken. The dependence of on A^ comes from the overlap S whose powers 
appear in the CCI equation (10), but, of course, are absent in the GP equation. As the 
number of bosons is increased, the shape of the CCI orbital is flattened and tends to the GP 
solution (which is constant at I7I = 0.2 for any A^ [7], see Fig. lA). Nonetheless, we see that 
the spike characterizes the CCI orbital for any number of particles, just as it does for the 
exact solution for N = 2, see Fig. 1. The same would hold true for other values of I7I and 
is a signature of the correlation energy accounted for by the CCI ansatz, and of the contact 
interaction U. 

In conclusion, within the CCI ansatz any eigenstate with a definite momentum of a 
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many-body Hamiltonian is exactly written in terms of a symmetry-broken function $. For 
two interacting fermions/bosons, $ has been determined exactly for all energy levels and any 
inter-particle interaction. For more particles and any given ansatz for $, the resulting CCI 
equation provides lower energies than the ansatz ^ = $ can. All of these makes our CCI 
ansatz a very attractive tool to be utilized in the many-body problem of an homogeneous, 
interacting fermion/bosons system, already with a Hartree-Fock/Hartree ansatz for <J>. We 
briefly remark that since the CCI ansatz is a variational one, it can be reformulated in the 
non-homogeneous case as well. 
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FIG. 1. The CCI orbital which exactly solves the problem of A^ = 2 attractive bosons on a ring 

For comparison, the corresponding GP orbital is presented: (A) I7I = 0.2; (B) I7I = 1.0. 
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FIG. 2. The CCI orbital for iV = 5, 25, 100, 1000 and 10000 attractive bosons on a ring for 
= 0.2. 
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